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Abstract 

^ In the context of survival analysis, [Marshall & OUdn] ( |1997| l introduced families of dis- 

^^ tributions by adding a scalar parameter to a given survival function, parameterized or not. In 

^^ that paper, we generalize their approach. We show how it is possible to add more than a sin- 

j£\ gle parameter to a given distribution. We then introduce very flexible families of distributions 

for which we calculate some moments. Notably, we give some tractable expressions of these 
moments when the given baseline distribution is Log-logistic. Finally, we demonstrate how 
to generate sample from these new families. 
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1 Introduction 



By various methods, new parameters can be introduced to expand families of probability distribu 
tions. This is an important issue in survival analysis ( |Lawless[|2003t Lee & Wang[ 2003 1 Marshall 



& OMn 2007 1. For instance, although the Weibull distribution is often described as flexible, its 



hazard function is restricted to being monotonically increasing or monotonically decreasing, or 
constant. Typically, in survival analysis, the limitations of standard distributions led naturally to 
interest in developing extended distributions by adding further parameters to a given distribution. 



The approach of Marshall & OMn ( 1997 1 works as follows. Let So{x) denotes the survival func- 



tion of a random variable, the corresponding Marshall- Olkin extended distribution has survival 

function, 

aSo{x) 



1 - (1 - a)So{x) 
where a > is an added scalar parameter. These distributions have been used in many areas: 



climatology: Biondi, Kozubowski, Panorska & Saito (2008 1; 



hydrology: Jose, Naik & Ristic (2008 1; 



insurance and finance: Garcia, Gomez-Deniz & Vazquez-Polo ( 2010| ), Jayakumar & Mathew 



(2008 1 and Kozubowski & Panorska (2008 1; 



medicine: Economou & Caroni ( 2007| ), Ghitany, Al-Awadhi & Alkhalfan ( 2007| l, Ghitany, 



Al-Hussami & Al-Jarallah] ( [2005) 1, |Gomez-Deniz| ( |2009] l, [Gupta & Peng| ( |2009l l and [Joser 



Naik & Ristic (20081; 



engineering: [Adamidis, Dimitrakopoulou & Loukas[ ( [2005[ ), [Adamidis & Loukas[ ( |1998[ l, 
[Gupta, Lvin, & Peng[ ( [2010| ), [Prabhakar Murthy, Bulmer & Eccleston[f2004[ ), [Silva, Barreto- 
[Souza & Cordeiro| ( [20T0[ ) and [Zhang & Xie[ ( |2007] l... 



For a review on the use of the extended Marshall and Olkin distributions, one can see [Nadara 



jah (2008 1. There is also a lot of methodological works around the Marshall and Olkin extended 



distributions. Independently of Marshall & OMn ( 1997 1, Adamidis & Loukas ( 1998 1 introduced 



the same type of extended distributions but only in the baseline WeibuU case. In Adamidis, Dimi 



trakopoulou & Loukas ( 2005 1, it is shown that the Marshall and Olkin extended distributions can 



be viewed as continuous mixtures. The basehne exponential case has been investigated by Silva, 



Barreto-Souza & Cordeiro (20101 and Srinivasa-Rao, Ghitany & Kantam (20091. For the more 



general baseline WeibuU case, one can see Ghitany, Al-Hussaini & Al-Jarallah (2005 1; Zhang & 



Xie (20071; Gupta, Lvin, & Peng (2010i. 



In Sankaran & Jayakumar| ( |2008| ), it is shown that the Marshall and Olkin extended distribu- 
tions satisfy the property of proportional odds functions. This property has been used notably by 
Economou & Carom] ( |2007[ ), |Caroni| ( |2008[ ) and [Gupta & Peng| ( |2009[ ) to introduce covariates. 



Some authors have investigated the behavior of the Marshall and Olkin extended distributions 
when the baseline distributions are Lomax, [Ghitany, Al-Awadhi & Alkhalfan ( 2007 1, Burr, and 



Pareto, Jayakumar & Mathew (2008), g- WeibuU, Jose, Naik & Ristic (2008 1 and normal Garcia, 



Gomez-Deniz & Vazquez-Polo (20101 distributions. Moreover, using the Marshall and Olkin ap- 



proach, G6mez-Deniz[ ( 2009] ) introduced a generalization of the discrete Geometric distribution. 



In the baseline exponential case, Kozubowski & Panorska ( 2008 ) highlighted the link between the 
truncated logistic distribution and the Marshall and Olkin extended one. 

Finally, some authors have investigated the generalization to the multivariate case, one can see for 
instance: [Thomas & Jose| ( |2004| ), [Sankaran & Jayakumar] ( |2008| ), [Jose, Ristic & Joseph] ( [2009] ), 



Yeh (2009 1 and Yeh (2010). Random minima and maxima related to Marshall-Olkin distributions 



have received attention in the litterature, |Amold| ( |1996[ ), [Biondi, Kozubowski, Panorska & Saito 
( |2008| ), [Ghitany, Al-Awadhi & Alkhalfaiij ( |2007l ), [Jose, Ristic & Joseph[ ( |2009l ), [Thomas & Jose 
( [2004| ), [Yeh| ( [2009] ) and [Yeh| ( [2010] ) . 





Figure 1: Hazard rate curve of a standard Weibull distribution with parameter (2, 2) (left), of an 
extended Marshall and OMn distribution with added parameter a = 0.05 (baseline Weibull (2, 2)) 
(center) and of an extended distribution obtained using our approach: two parameters (oi = 10~^ 
and a2 = 0.15) are added to a baseline Weibull (2, 2) (right). 

In that paper, we generalize the approach of Marshall and OUdn in order to obtain more flexible 
distributions. Indeed, coming back to baseline Weibull one, it is easy to show that the hazard rate 
curve of the extended Marshall and OMn distribution has at most two waves. As shown in Figure 
[T] using our approach, we can obtain more rich hazard rate curve. Our proposal is introduced in 
Section [2] In Section [3} we investigate some properties of these extended distributions for various 
baseline cases. 

2 The new semi-parametric family 



In Marshall & Olkin (19971, it is introduced a new semi-parametric family of probability dis- 
tributions. One real parameter a > is added to a given univariate probability distribution 
Fq : M — ;• [0, 1] by using an increasing function denoted by Qa- Indeed, the Marshal and OMn 
extended probability distribution is given by 



9a{Foix)) 



Fo{x) 



a + {l-a)Fo{x) ' 



Here, we generalize this support function with an arbitrary number of external parameters. 



Let 5o = 1 - -Fo and q G N*. Clearly, for any (ai, . . . , a^) G ]0, +oo['^ and any u G [0, 1], we 
have, 

We define the function gai,...,aq ■ [0, 1] — )■ M by 



q 

u 

i=2 



Y[{ai + u- aiu) 



9ai,...,aq{u) — q ^ q / q \ "l '^ ' 

As a generalization of the Marshal and OMn proposal, we define the Let us now define the func- 
tions -Fai,...,ag and Sa^^...^aq such that : 

Fai,...,aq : M — > M 

Theorem 2.1 For any integer (7 G N* and any (ai, . . . , a^) G ]0, +oo['^, the application Fa^^,,,^aq 
from M to M takes its values in [0, 1], and the function Fa^^,„^aqfrom M to [0, 1] defines a probability 
distribution. 

Proof 

The proof is given in Appendix. D 



If we replace in Theorem 2.1 the cumulative probability function Fq by the survival function 5*0, 
we still obtain a survival function. This defines a new probability distribution. In that work, we 
only consider the initial construction. 



For any a > 0, the function Fa,..., a from M to [0, 1] is a Marshall-Olkin extended probability 
distribution. That is, for any real number x, we have, 



or, in another way. 



^a,...,a\^) 



a + {l-a)Fo{x) 

aSo{x) 
1 - (1 - a)So{x) 



Let us now introduce the following lemma. 

Lemma 2.1 We suppose that there exists t G R such that < Fo^t) < 1. 

In that case, there exists a unique interval JofM such that, for any x £ (R — J)n] — oo, t] we 
have Fo{x) = 0, for any x G (R — J) n [t, +oo[ we have Fo{x) = 1, and for any x in J we have 
< Fo{x) < 1. 

The proof of this lemma is trivial. The interval J does not correspond to the usual definition of the 
support of a probability distribution (the smallest closed set whose complement has probability 
zero). We use that specific definition for technical reasons. 



2.1 



Let us suppose that Fq is continuous on R and C^ on J, the interval J is defined by Lemma 
Let /o be the probability density function of Fq. We can take /o as, 

/o: R ^ R+ 

X I-)- Fq{x) if X G J 
otherwise 

Here, it is very easy to verify that for any q £ N* and any (oi, . . . , a^) G ]0, +oo['', the function 
Fai,...,a„ from R into [0, 1] is continuous on R and C^ on J. 



Moreover, the probability density function fai,...,aq of Faj,...^aq can be taken as, 



Jai,...,aq ■ 



X ^ Fa^,...,aqi.x) = 9au...,aq {Mx)} fo{x) if X G J 

otherwise 



For instance if g = 2, we get for any x £ J, 



_ ^ aia2Fo{x) + aJFoix) - 0102 - aj - 2a2Fo{x) ^ 
{aiFo(a:) + a2Fo(x) - ai - 02 - 2Fo{x)} 



Let us now introduce two Lemmas. 

Lemma 2.2 Let g G N* and ai, . . . ,aq be strictly positive real numbers such that 
ai + ■ ■ ■ + aq > |. For any n G [—1, 1], 

9 / <? \ 



i=i \i=i 



and 



where, for any m G N*, 



9ai,...,aq{u) — 2_^ Cq^ai,...,aq 



m=l 



Q' nL2«* 



•'q,ai,...,aq,m 



with 



('^-i)KELi«on,,,.f:l,,,, ' 



1 






CTi = < 



if i = 



^2<ji<...<ji<q aj^ 




"^ X • • • X — ^ if 1 < i < o - 1 



if i > 



Note that we have 

oo 

/ J (^q,ai,...,aq,m = -L 
m=l 

Proof 

The proof is given in Appendix. D 



We trivially deduce the following corollary used in the proof of Theorem 2.2 



Corollary 2.1 Let g G N* and ai, . . . ,aqbe strictly positive real numbers such that ai+- ■ ■+aq > 
q and Oj < 1/or 2 < i < q. Then, for any m G N*, 



Cq,ai,...,aq,m ^ U • 



Lemma 2.3 Let g € N* and ai, . . . ,aq be strictly positive real numbers such that 

ai + • • • + Qq < 2q. For any u G [0, 2], 

^ffli- ["^ai- q\u^O, 



i=i \i=i 



and 



where, form G N*, 



gau...,aqiu) = 1 + X] dq,a^,-,aq,m{u - 1) 



m=l 



"9 n-^m 



d,,a„...,a„m= ^ {{m + q - 1) ■ ■ ■ {m + l){^^f^] 



with 



(Ti = < 



Y.2<h<...<n<n{^ "-jl) ^ 



if i = 
••x(l-a,J if l<i<q-l ' 







if i > 



Proof 

The proof is given in Appendix. D 

We now study the expectation of our family of distributions when the baseline probability distri- 
bution Fq is continuous on M and C^ on J. 

Theorem 2.2 Let {^, A, P) be a probability space, q eW, {ai,...,aq) g]0,+cx)['^, 

Xq : il — 7- M and ^ai,...,ag : 0, ^ M. be random variables such that Xq ~ Fq and Xa^^^^^^aq ~ 

Fai,...,aq and w :M ^Mbe a borelian function. 

If Fq is continuous on M and C^ on J and if E{\w o Xq\) = J^ \w{x)\fQ{x)dx G M4-, then 



Moreover, ifai,...,aqbe strictly positive real numbers such that ai + ■ ■ ■ + aq > q and Oj < 1 
for any 2 < i < q, then 

E{\w o Xa,,...,aJ) < aiE{\w o Xo\) . 

Proof 

The proof is given in Appendix. D 



As a direct consequence, if ai, . . . , a^ be strictly positive real numbers such that ai + • • • + Og > g 
and a-i < I for any 2 < i < q and, f^ \x\^ fo{x)dx G M+ for a real number r, then 

E{\Xa„...,aS)<a,E{\Xon. 

The following proposition will be used to derive explicit formulae on some expectations when the 
baseline probability distribution is a classical Log-logistic. 

Proposition 2.1 We suppose that Fq is bijective from J to ]0, 1[. Let (17, A, P) be a probability 
space, (ai,a2) € ]0,+oo[ such that oi / 02, -^ai,a2 : i^ — > K and V"ai,a2 '■ ^ ^ ^ be 

Define for any w G il 5'mc/i that Xa^^aii^) ^ <^> 



|ai-a2|r 



^ai,a2(w) =log 



1 



Fo{^ai,a2(^)} 



1 



and for any lj £ Q such that ^ai,a2 (i^) G IR — ^, 



y,,,a,H = o. 



Then the random variable, 



^ai,a 



^ai,a2 + Ki,a2 " log ^^ if Oi > 02 



^ai,a2 I ''ai,a2 -^OJ 



ai+a2 



if ai < 02 



/i distributed according to a logistic distribution. 
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Proof 

Let us denote by aai,a2 the cumulative probability function of T4i,a2» by r^^ ^j the one of l'ai,a2' 
by ta^^a2 the one of -Fai.aa- and by Ua^^ai the one of Lai,a2- 
For any u G R, we have 



^ai,a2V^j — J^ 9ai,a2 



1 

1 + exp(u) 



and 

exp(ii) 



tai,a2\'^) — 9ai, 



'"^ ' 1 + exp(ii) 
We first suppose that ai > a2 and for any {u, t) G R^, we define 



2(ai+a2)exp(^+^t) 

ii,a2,w(,*j 



{2 exp(t) + (ai + 02) exp(ti)} 
For any {u,t) G M x [0, +00 [, we clearly have 

6'ai,a2,«(*) = V,a2(^ " *)"ai,a2(*) 

and 

/ / {u-t)a' (t)dt= 2(«i+«2)exp(n) 

/ 'ai,a2l" ''J«ai,a2V'';"'' ro , ^ , \ / ^^2 • 

JR+ |2 + (ai + 02) exp(u) j 

Therefore, for any w G M, 

/ / r' a{u - t)a' {t)dtdu = -—j — -^— — = lia,,^^ ( w - log — - — ) 

J-ooJr+ 2 + (ai + 02) exp(t;) V 01 + 02/ 

and 

exp(v) 
1 + exp(f ) 
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We now suppose that 02 > ai and for any (n, t) G M? and we define 

2(ai + aa) exp(u - -^t) 



122— 11 



{2 exp(u) + (ai + 02) exp(f )} 
For any {u,t) G R x [0, +00 [, we clearly have 

6'ai,a2,«(*) = *ai,a2(^^ " *)«ai,a2(*) 

and 

,, , ,. , ,,w, 2(ai +a2)exp(M) 

*ai,a2(" - 0«ai,a2(*)'^* = J^ T^- —^2 ' 

{2exp{u) + ai + 02) 
Therefore, for any f G M, 

" / ,/ / ,x / .,N,,, 2exp(w) / ai+a2 

ooyR+ 2exp(t!) + ai +a2 V 2 

and 

exp(f) 



Uai,a2{v) 



1 + exp(f ) 

D 



In Marshall & Olkin ( 1997 1, it is shown a specific property of the introduced semi-parametric fam- 
ily: the proposed distributions are geometric extreme stable. That is, they are both minimum and 
maximum stable, when the random indexation variable of the sample size is distributed according 
to a geometric distribution. We now consider maximum stability for our semi-parametric family. 



12 



Proposition 2.2 Let g G N* and ai, . . . ,aq be strictly positive real numbers such that ai + • • • + 
o-q ^ Q cmd ai < Ifor any 2 < i < q. Let {0,, A, P) be a probability space and Na-^^^^^^aq '■ ^ — ^ 
N* be a random variable such that, for any m £ N*, 

Pr(iVai,...,aq =m) = Cg^ai,...,aq,m ■ 

• If, for i G N*, Xi : VL ^ M. be mutually independent random variables, independent from 

Nai,...,aq, such that Xi ~ Fq, then 

Ki,...,a, = max(Xi, . . . ,X7V„^ _„J ~ Fa^,...,a„ ■ 

• If, for 6 > and i G N*, Xi : Q ^ M be mutually independent random variables, indepen- 
dent from iVaj^...^a, such that Xi ~ -Fb,...,b, then 

Va^,...,aq = max(Xi, . . . ,X7V„^ _„J ~ i^feai,...,fea, • 

Proof 

Part one: we have for any real number x. 



pr(V;i,...,a, < X) = ^ pr(iVai,...,a, = m-)^o( 



X) 



m=l 



pr(Vai,...,a, <X)= ga^,...,aq {Fo{x)} , 
Pr(Va,,...,a, <X)= Fa^^...,aq{x) ■ 



13 



Part two: for any (u, ai, . . . , a^, b) G [0, 1] x ]0, +oo[''''' , we have, 

9ai,...,aq {gb,...,biu)} = gbai,...,baqi^) ■ 

D 

More generally for q >2, the new probability distributions are not maximum stable. 

3 Some properties of the extended distributions 

3.1 Fo is WeibuU 

We first suppose that Fq{x) = (1 — exp{—x))Ix>o (exponential distribution with expectation 
equal to 1). 

Lemma 3.1 Let Fq{x) = (1 — exp(— x)) Ia;>o <^f^d /o is the corresponding probability density 
function. For any r > and any m G N*, we have 

l^x^Foixr~'foix)dx = rT{r)Y, (^ J ^)(-iy(i + 1)"''"' • 



14 



Proof 

Let r > 0, p G N and z > 0, 



and 



Then, 



and 



Fo(zf = X:(^)(-l)^exp( 



-jz) 



z^FoizfMz) = ^ (P\-iyz-exp{-{j + l)z} . 

Fo{xrfo{x)dx = y2(^)i-^y [ x^-expi-a + iM^x 



/ x^Fo{xrfo{x)dx = J2 (^) (-i) \.'\^'L = ^r(o X] (^) i-iyu + 1 



-r-1 



D 



Proposition 3.1 Let Fo{x) = (1 - exp(-x)) Ia;>o and Xa-,,...,aq ~ -?^ai,..,a,- 

//'ai , . . . ,aqbe strictly positive real numbers such that ai + • • • + a^ > q and ai < I for 2 < i < q, 

oo m—l / ^\ 

E{Xa,,...,a:) = rT{r) ^ mc,,,„...,,„^ T. ( ~ ) (-l)'(-?' + l)"''"Vor an); r > 0. 
Ifai,...,aqbe strictly positive real numbers such that q < ai + ■ ■ ■ + aq < 2q and a^ < Ifor 

oo 

2<i<q, E{Xa,,...,aJ) = rT{r) ^ {-ir'Uq,a„...,a„mm-^ for any r > 0. 



m=l 



Proof 



The proof of this result is almost trivial using Lemmas 2.2 2.3 and 3.1 D 
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We now consider the case of the classical Weibull distribution with two parameters, that is 



Fo{x) = F;;"'^x) = l-exp 



b2' 



"■x>0 



^bi,b2 



j^biM / 7761 ,62 



with (61,52) G]0,+oop. Let Xa\\.'^.,aq ~ Fall'.^.,aq {Fall^,ag bc the corresponding parameter 



biM 



'1,1 



augmented distribution). We can easily verify that Xal]. .,a„ = biXa[,...,a„ in distribution and, 



then for any r > 0, 



^ai,...,aq 



i^te''^/) =b,'E Xl'} 



Finally, we consider the case of the generalized Weibull distribution with three parameters, that is 



^6l,fe2,b3(^) 



1 — exp 



, X b2^ b-i 



llx>0 



3 T ^ vbi,b2,b3 _ 7761,62,63 

Xal','!^.',al = bil{l + Xa[\'.t,aqy'' — 1 ^ ^ in distribution and then, if we suppose in addition that 
7712 = r" G N — {0} and 63 G N*, it follows that for any integer m G N*, 



with (61, 62, 63) G]0, +oo[3. Let X'all'^fal ~ F^ll'X- We can easily verify that 

-61,62,63 _ J, f^i , x-1,1,1 \b-i ^ 1 6^ 



E{x^'!:^^^ 



mm2 b:jk 
fc=0 j=0 



(7^)('f)(-l)"— 'i5(xi;;:',„;*") 
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3.2 Fq is Log-logistic 

We now consider the case of the standard Log-logistic distribution, that is 



^o(x)=(^)I.>o^ 



Proposition 3.2 Let Fq{x) = [rhj ^x>o and Xai^...,ag ~ i^ai,...,a,- 

Ifai, . . . , Gq be strictly positive real numbers such that ai + - ■ ■ + aq > qandai < lfor2 < i < q, 
then for any r such that \r\ < 1, 

oo 

-E'(^ai,...,a/) = ^ "iCq,ai,...,ag,mBeta(l -r,m + r) . 

m=l 

Ifai,...,aqbe strictly positive real numbers such that q < ai + ■ ■ ■ + Uq < 2q and ai < Ifor 
2 <i < q, then for any r such that \r\ < 1, 

CX) 



m=l 



Proof 



The proof of this result is almost trivial using Lemmas 2.2 and 2.3 □ 



We now consider the case of the classical Log-logistic distribution with two parameters, that is 

with (61,62) G]0,+oop. 
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,b2 



Let Xal].^.,a„ ~ ^ai,'..^,a„- As in the WeibuU case, we can easily to verify that Xa\].^.,a„ 



biXa[,...,aq '"^ in distribution, and then for any r such that \r\ < 62, 



^(^a';:'.%;) = b/Eixii' 



Proposition 3.3 Let Fq f j^ j Ia;>o and Xa^^a2 ~ -^01,02- For any r such that \r\ < 1, 

EiXaua2 = 5 ■ . X r + 1 . 

\ 2 J sm(r7r) \ai + a2 J 

Proof 

When the baseline distribution is a classical Log- logistic with two parameters, the extended Marshall- 
Olkin distribution with one external parameter gives again a classical Log-logistic distribution. 
Therefore, the case oi = 02 is trivial. 

Let us now consider the case ai > 02 



Using Proposition 2. 1 we have 

ai + a2 



exp{-Laj^^a2) = exp(-y;^,a2)exp(-V;^ 



(12 J 



Then 



and 



' {exp{-La^^a2)V = {ew{-Yai,a2)V {e:^P{-Vai,a2)y 



"^ + °' ' i^({exp(- V,,,)}^) = Ei{expi-Y,,,a,)y)E{{expi-Va„a2)r) 
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The random variable exp(— L^^ (12) has a Log-logistic distribution and for any r such that \r\ < 1, 



E{{expi-L,,,,,)Y' 



sm rvr 



Moreover, 






/o «! - «2 (oi - a2)r + (ai + 02) ' 

Therefore, 

'ai + a2V rvr (ai - a2)r + (ai + 02) 



Ei{eM-Ya,,a,)V) = EiXa„a/) = [^^^ 



2 J sin(r7r) ai + 02 

We use the same type of reasoning for the case ai < 02- □ 

If Xa\]a2 ~ fallal then for any r such that \r\ < &2 

3.3 Random sample generation 

When the distribution Fq admits a density with respect to Lebesgue measure on R, if we are able 
to generate a sample from Fq, we can use the accept-reject algorithm to generate a sample from 

^ ai,...,aq- 
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Indeed, following the proof of Theorem 2.2 for any (ai, . . . , Og) G ]0, +oo['', 



Jai,...,ag l^J ^ I min(l,ELi «>/9)' ' ^*=1 * ^' mm{l,ELi "./?)"+' 

Moreover, if ai + . . . + a^ > g and if for any integer i such that 2 < i < g, we have ai < 1, then 

/ai,...,a,(a;) < ai/o(x). 

4 Conclusion 



We shown how to generalize the approach of Marshall & OUcin (1997 1 in order to obtain more 



flexible families. We investigated some properties of the introduced distributions. We are now 
working on the parameter estimation task. The preliminary results are extremely encouraging, 
notably the application of our proposal on some real survival datasets. 
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Appendix 

Proof of Theorem 12.11 

The case q = lis trivial, it corresponds to the Marshall and Olkin derivation. 

Part one: 

Let us first prove that for any (oi, 02) G ]0, +oo[ , the function gai,a2 from [0, 1] into R is increas- 
ing on [0, 1] and takes values in [0, 1]. 

For any u G [0, 1], we have, 

4{aia2U + afu — aia2 — af — 2a2u) 



a'auaii^) 



(aiu + a2U — ai " a2 — 2u)^ 



Moreover, 



and. 



Therefore, 



aia2{u — 1) + a'l{u — 1) — 2a2U < , 
ai{u — 1) + a2{u — 1) — 2n < , 

{01 (n - 1) + a2{u - 1) - 2m}^ < . 



Thus, the function ga-^,a2 from [0, 1] to M is increasing on [0, 1]. As gai,a2 (0) = and gai,a2 (1) 
1, the function gai,a2 from [0, 1] into M takes values in [0, 1]. 
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Part two: 

Let us now prove by induction that for any integer q such that q > 2 and any (ai, . . . , Og) G 
]0, +00 ['^, the function gai,...,aq from [0, 1] to R is increasing on [0, 1] and takes values in [0, 1]. 

For any integer q such that q > 1, we denote by (Hq) the condition such that for any (ai , . . . , Og) G 

]0, +oo['', the function gai,...,aq from [0, 1] to R is increasing on [0, 1] and takes values in [0, 1]. 

Let q be an integer such that q > 3 and suppose that the condition (Hq^i) is verified. 
In that case, for any {u,ai,. . . , aq) G [0, 1] x ]0, +oo['^, we have, 

gai,...,aq{u) = g<n ^2^ 1 {gag,...,aq{u)} 



and. 






ga^ ^azii 1 \u) 



E£if: + i-{E£i^-(9-i)} 



1 ' 



(9-1) 



— Y5-l "9-1 l^j 

ag '*■■' ag 



E£if:-{E£i^-(9-i)} 



9-1 



E£i^ + i-{Eni^-('?-i)} 



u 



9 • 



Also, 



■q-l (q-l -\ 



=1 "9 



,1=1 ^ 



E 



■9-1 



-7-1 



E?- E 



,t=i ^ 



(<?-!) 



ti 



, i=l 



q-j 



Mj-l a, 






^<?-l flj 



(g-1) U 



^y-l Oj 



i=l a„ 1 ^ 



■ii=l a„ 



(g-1)} 



U 



<? 

9-1 ^ 2^ 
j=0 



■9-1 



■q-l 



E!^- E 



• 1 Og 



1=1 ^ 



{q-l)}u 



i-i 
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For any u G [0, 1], let us define, 



rjau...,a,iu) 



E£i^ + i-{E£i^ -(</-!)} 



u 



E£i^-{E£i^-('/-i)} 



u 



q-1 



We have, 



Va^,...,aS^) 
-9-1 



|s-<-+§C)<-^>{fe-"->}[fe-{fe-'' 



i)U 



'9-1 



E^-(^-i) i + E (1-^- 



.1=1 



i=2 



■g-l 



-g-1 



n ~J 



Ez- E^-(^-i 



1 = 1 ^ 



.4=1 ^ 



and then. 



^ai,...,a,(l) 



{|^<.-.}{'^|:(;')('-K,-:) 



For any a; G R — {0}, let us define. 



i=2 ^-^^ 



-i+i 



We have 



rq(x) = X 



+ 1 -1 



rg(x) = X - + 1 



1 



X ■ 
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and, 



rg(x) 



i=2 ^-^^ 



-^i + iV-^fi + i 

X J X \x 



9-1 



+ 1 



9-1 



l-q 



+ 1 -1 



Therefore, 



Moreover, we have. 



r'Jq-l) + l = 0. 



<,...,a,(l)=0. 



Vai,. 



■q-1 



(9-1) 



E 

i=2 



(i-i)i 



■g-i 

E 

• 1 ""a 
.1=1 ^ 



an 



■q-l 

E? 



{q-l) 



u 



-i-1 



Then clearly, <,...,. (n) > and r/^^^ 



— ^ai a (1) = 0- Therefore, the function 



li I— ;• 1 — r^ from [0, ll to R^ is increasing on fO, 11. 

Using condition (Hg^i), the function tt i— )■ (701 a^-i (u) from [0, 1] to M is increasing and posi- 
tive. It follows that the function gai,...,aq from [0, 1] to M+ is increasing on [0, 1]. 
Moreover, gai,...,aq{0) = and 5'ai,...,a, (1) = 1- Thus, the condition (Hg) is verified. 
Therefore, we get the result by induction. 

Part three: 



Let us now prove that the function Fai,...,aq from M. to [0, 1] defines a probability distribution. 
The function Fa-^^...^aq from M to R takes its values in [0, 1], and it is a right continuous increasing 
function. Moreover, lim^^_oo Fa^,...,aq{x) = and lim^^+00 Fa-^,...,aq{x) = 1. 
Thus, it is a probability distribution. 
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Proof of Lemma 1121 

Let (ai, . . . , Og) G ]0, +oo['^ such that ai + • • • + Og > |. For any u G [—1, 1], we have, 






Then, 






Moreover, 



«fl{«. + (1 - «^M = (fl«.)^fl(l + ^-^^n) = {\\ai)u\{{l + ^-^^u) , 

i=2 4=2 i=2 j=2 i=l ^ 

g g g-1 g g-1 g g g oo 

-« JJ {fflj + (1 - ai)u} = ( JJ aj)^ ^ aiu' = ( JJ Oj) ^ aiu'^^ = (JJ Cj) ^ o-j_in* = (JJ Oj) ^ aj-iu^ 

1=2 i=2 4=0 4=2 4=0 4=2 4=1 4=2 j=l 

Therefore, 

W ^J- 12^4=1 «*J ,n=l I j+fc=m,l<i<g ^ 2^4=1 "« / 
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Proof of Lemma 12.31 



Let (oi, . . . , Cq) G ]0, +oo['^ such that ai + ■ ■ ■ + Uq < 2q and u G [0, 2]. We define w = u — 1, 
we get that w G [—1,1] and, 



yai,...,aq\,u,) — q 



9 ' 



{ELi«i-(ELi ««-«)"} 



{T,i=iai-{Ei=iai-q)il+w)y 



5ai,...,a,W 



(i+"')nL2{i + (i-«*M 






(1 - "'=t ' ' w 



Moreover, -1 < ^'-i"' "^ < 1, then -1 < S=L^^^y; < i and, 



(1-^H^^.r ('^-D't^o 



(fc 



+,-i)---(^+i){^^^i^^^y 



Also, 



g-l oo 

m,=0 r7i=0 



n {1 + (1 - «.>} = E ^-^" = E 

i=2 m,=0 r7i=0 

q oo 

(1 + tu) J]^ {1 + (1 - aO^i^} = Y, ^'"^'^ + Z] ^^ 

i=2 m,=0 m=0 

? oo 

(1 + u;) J]^ {1 + (1 - ai)u'} = 1 + ^ (a^ + cj„_i)u; 



oo 



m+1 



1=2 



m=l 



Then, 



gau...,aAu) 



(q 



rw 1 1 + X] (^"^ + ^'»-i)'(i'™ > 



oo 

E 

.fc=0 



(fc+,-l)---(^+l){^^-i^^^}' 
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9au...,agiu) 



(<? 



ky. {Er=o(^ + q-l)---ik + l)(^t^) V 



+ Em=l Efc+i=m,l<i<g(^ + ^ - 1) . . . (/c + 1)( 



ELi^'-g^ 



(o-j +(Tj_i)t(; 



5ai,...,a5('UJ 



1 + 



(5-1)! 



E::=l{(m + <?-l)...(m+l)(^^^: 



ELi°'-g ^ 



+ Efc+,=„.,i<,<,(fc + q-l)...{k + l)(^-f^) (CT, + a,_i) 



w 



Proof of Theorem O 



Part one: 



Let us prove first that for any integer q such that q > I, any (ai, 
u G [0, 1], we have, 

5ai,...,a„(l) =«i> 



,aq) G ]0,+oo['' and any 



ffll,...,a,(^)>0: 



< 



ni'=2 niax(l,a, 
min( 






Oi-gl 



5ai,...,a,(«) 

n-=2 max(l, ai) EL2 |1 - «il n2<j<gjVi max(l, 



min(l, ^'=^ ' ) 



+ 



min(l, 



ELi °' ■ 



For any u e [0, 1], we have, 



5ai,...,a,(w) 



{ELi ai-(ELi a.i-<i)uy 

"ry lZ_^j=l "« 'i'l^q ^^^^9 N 19+1 

{ELi ai-(ELi Ci-g)" j 
{ELi o^i-CELi '^'-'j)"}'' 



27 



Then, 

9'au...,a,i^) = 1 + J]ai-g + J](l-ai) = ai. 

i=l i=2 

For any u G [0, 1], we have 

9ai,...,a,W^ 19 {ELi«.-(ELi«.-gH'' 

{ELi "»-(ELi "«-'?)"}'' 
Moreover, for any integer 2 < i < q 

< min(l, tti) < tti + u — atu < max(l, ai) , 

q q q q 

< min((7, N^ Oj) < >^ Oj — ( >^ Oj — g)u < max(g, N^ Oj) 



i=l i=l i=l 1=1 



Therefore, 



n^=2™'^^(^'"') 



/ / N ^ n \\- oniax(l,a7 



,„g+l|V« „. „|Jli=2iWMi) 
' -^ mm(<j,X,i=i Qi) 

^q EL2ll~°'in2<j<9.J7^ir"ax(l,aj-) 

5 ;rvj ) 






yai,...,a„l"J ^ Y-9 „. g"r|/,"« yl ^q a+l ~'~ Y"« « 5 

mm(l,^^) tr min(l,%^)' min(l, ^^) 

As the function gai,...,aq from [0, 1] to M is increasing and differentiable on [0, 1], for any u G [0, 1], 



5ai,...,a,(«) >0. 
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Part two: 

Let us suppose that f^ \w{x)\fo{x)dx G M+. We have 

E{\w o Xa,,...,a,\) = / \w{x)\g'^ {Fo{x)} fo{x)dx . 

Using the resuh given in Part one, for any u £ [0, 1], 



[min(l, '-1 ' ) 



+ 1 2^i=l ^i ~ 1\ ^ ^ g+T 

min(l,5i^) 

min(l, '=1°' ) J 



Therefore, 

^(ko^ai,...,aJ)GM+. 

By the corollary |2.1[ if g > 1 and ai , . . . , Og be strictly positive real numbers such that ai + • • • + 
ttq > q and Oj < 1 for any 2 < i < q, then we have that for any integer m which is greater than or 
equal to one and any u G [0, 1], 

<,...,a,(^)>0 



Therefore, 



Sa^,...,ag{u) < 5ai,...,a,(l)i 

E{\w o Xa,,...,aJ) < aiE{\w O Xo 
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